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The event-by-event azimuthal fluctuations in high-energy heavy-ion collisions are analyzed by 
means of the so-called measure. The fluctuations due to the collective transverse flow and those 
caused by the quantum statistics and resonance decays in the hadron gas are discussed in detail. 
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Various phenomena can lead to the nontrivial azimuthal fluctuations in (ultra-)relativistic heavy-ion collisions. One 
Q-T mentions here jets and minijets being the result of (semi-)hard parton-parton scattering [|lj and the collective transverse 
flow due to the anisotropic pressure gradient ||||- The latter phenomenon, which is well known at the intermediate 
collision energies (for a review see e.g. Q), has been recently observed in the nucleus- nucleus collisions at CERN SPS 
§,§. As argued in the series of our papers the color plasma instabilities also generate a transverse collective 
flow in heavy-ion collisions at RHIC and LHC. Since the final state azimuthal fluctuations appear as remnants of the 
inhomogeneous or anisotropic early stage of the collision, when the minijets are copiously produced (|] or the color 
instabilities occur Q, the effects of interests are expected to be rather small ||. Thus, it is a real challenge to extract 
Q\ . them from the statistical noise. 

A method, which has been successfully applied to study the collective flow in heavy-ion collisions at SPS energies 
H , is based on the Fourier analysis of the azimuthal distribution [l^] . Specifically, one expands the distribution 
of the particle azimuthal angle measured with respect to the reaction plane, which is reconstructed on the event-by- 
. event basis, into the Fourier series. Then, the first harmonic tells us about the so-called directed flow, the second one 
' about the elliptic flow, etc. Although the method is very powerful there are subtleties in applying it to the real data. 
In particular, the finite statistics and non-flow interparticle correlations - the jets are expected to be an important 
,—H ' source of such correlations at RHIC - distort the reconstructed event plane [|llj. Due to the finite resolution of the 
reconstruction procedure, it is difficult to measure the Fourier coefficients higher than the second one |p^| . 

The Fourier expansion method has been designed to study the collective flow which is correlated with the reaction 
plane. However, the jets or color instabilities, which have been mentioned above, generate the transverse collective 
motion being, at least approximately, independent of the reaction plane orientation. Then, the Fourier expansion 
method, as developed in J5~0| — |T"^] , is not sensitive to such phenomena and one needs another tools for the data analysis, 
^ ' for example the method which focuses on the two-particle large-angle azimuthal correlations |l3| . 

In this letter we propose to study the azimuthal fluctuations by means of the previously introduced fluctuation 
measure $ |Q. Although the measure has been invented for other purposes, see below, it might be also useful for 
the case of azimuthal fluctuations. The measure analysis does not demand the reaction plane reconstruction and 
can be easily applied to the data. <E> is also sensitive to the different sources of correlations. It appears difficult to 
disentangle various contributions but the integrated information provided by <E> can be combined with that offered 
by other methods, in particular the Fourier expansion. For example, we show in this note that the measure of 
(^—fluctuations deviates from zero when the transverse collective flow is present in the system. Since all harmonics 
contribute to $ one can check whether the first two Fourier coefficients, which are measured, saturate the observed 
value of $. The analysis of measure of (^—fluctuations can be also combined with the $— measure studies of 
other kinematical variables. of pr~ fluctuations has been already shown to be very sensitive to the presence of the 
jets Thus, the simultaneous measurement of 0— and pt~ fluctuations can be helpful in disentangling various 

contributions. One should remember however that except the dynamical correlations caused by the hydrodynamic 
flow, jets or color instabilities, there are also correlations due to the quantum statistics of particles. These correlations 
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are always present. Therefore, we compute $ in the ideal pion gas to estimate how the Bose statistics influences the 
azimuthal fluctuations. We also consider here the effect of hadron resonances. 

Let us introduce the correlation (or fluctuation) measure <&. We define the variable z = x — x, where a; is a 
single particle characteristics such as the particle transverse momentum or the azimuthal angle. The overline denotes 
averaging over a single particle inclusive distribution. In our further considerations, x is identified with the particle 

azimuthal angle. The event variable Z, which is a multiparticle analog of z, is defined as Z = Yl^L^ {xj — x), where 
the summation runs over particles from a given event. By construction, (Z) = 0, where (...) represents averaging over 
events (collisions) . The measure $ is finally defined in the following way 

(1) 

As shown in (TJ], $— measure possesses properties which are very useful in the data analysis. It equals zero when 
the final sate particles are independent from each other i.e. there is no inter-particle correlation in the system. The 
second property is less trivial: the value of <!> is exactly the same for the nucleon-nucleon (N-N) and nucleus-nucleus 
(A-A) collisions, if the A-A case is a superposition of N-N interactions with no secondary scatterings. In other words, 
the strength of the correlation is not influenced by the number of uncorrelated particle sources. The measure $ has 
been successfully applied to the NA49 experimental data and it has been found that the dynamical transverse 

momentum correlations, which are present in N-N collisions, are significantly reduced in the central Pb-Pb reactions. 
The correlations observed in these collisions are fully explained by the effect of Bose statistics of pions |l8],|o| . 

Let us compute the $— measure of the azimthal fluctuations caused be the transverse collective flow. The inclusive 
4>— distrubtion is, of course, flat i.e. 

Pinc{4>) = ^ 6(0) 6(2tt - <f>) , (2) 

which gives </> = ir and </> 2 = §7r 2 , and consequently, z 1 — ■ Following |icj , the azimuthal distribution of a single 
event is chosen in the form of the Fourier series i.e. 

Pev((f>) = y I 1 + Z! ( x nCOs(ncj)) + y„sin(ra0))] 0(0) 6(2vr - 0) , (3) 

71=1 

where the parameters X n and Y n change from event to event. The distribution (^) is usually rewritten as 

P«,{4>) = ^- [l + 2^u n cos(n(0-^ n ))]e(^)e(27r-^) , (4) 

n=l 

where 



X n = 2v n cos{mp n ) , Y n = 2v n sm(nip n ) ■ 

The first term of the Fourier series from eq. (H) or (|]) corresponds to the so-called transverse directed flow. The 
amplitude vi controls strength of the flow while %p\ determines the reaction plane. The meaning of v n and ip n for 
n > 1 is analogous. 

The angles tp n vary from event to event and their distributions are flat. Therefore, the event distribution (Q) averaged 
over events provides, as it should, the inclusive distribution (0). Further, we consider two extreme cases. In the first 
one, which seems to be appropriate for the hydrodynamic flow analysis, the ip n angles are maximally correlated to 
each other and uniquely determined by the reaction plane angle ip r i.e. ip r = nip n + ot n . Then, the averaging over 
events corresponds to the integration over the angle ip r - In the second case, the angles ip n are independent from each 
other and one integrates over all ip n to average over events. 

Since Z = ^2f =1 (<j>i — </>) one finds in the first case that 

<Z 2 )=— / dl/j r J2' PN d ^ / d<t>N Pev{4>l) Pev{<t>N) (01 + + 4>N ~ Nj) 2 , (5) 

where Vn is the multiplicity distribution. The formula analogous to (||), which corresponds to the second case, 
includes the averaging over all angles ip n . It is understood that there is one more averaging in eq. (p|) which is not 
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explicitly shown. Namely, the averaging over the amplitudes v n which also change from event to event. To simplify 
the notation we also neglect here a correlation between the event multiplicity and the flow strength. After elementary 
calculation, one finds from eq. (||) that 



with 



and 



(Z 2 ) = Y (N) + {(N*)-(N))S, 



N 



x ■ iXT=i (lf) 2 > for l-stcase, 
(( EZi V ir) 2 ) fOT 2-nd case. 



Finally, we get 




As expected, $ = for S = 0. When S — ► 0, (A^ 2 ) = (N) 2 and (N) ^> 1, we get an approximate expression 



<I» = JL {N) S . (6) 



Let us consider how large is the expected effect. The amplitudes v\ and v-i observed in Pb-Pb collision at 158 GeV 
per nucleon do not exceed for pions the value of 0.03 We take v\ = V2 = 0.03 and v n — for n > 2. We also 
neglect here the variation of v\ and V2- Then, one finds from Eq. (^J) that for (N) = 170 Q $ equals 0.058 in the first 
case and 0.105 in the second one. 

As already mentioned, the transverse flow is far not the only source of the azimuthal fluctuations. In particular, 
the quantum correlations contribute here. We compute $ in the ideal quantum gas to estimate the effect of quantum 
statistics. Modifying our previous calculations [p8|-p0[, one immediately finds 



(Z 2 ) 1 f d3 P fJ . -2 A"V E 



(8) 



(N) pj (2tt) 3 vr YJ {X- 1 ef JE ± l) 2 ' 
where 

d 3 p 1 



(2tt) 3 \- 1 ef" s ±l 



(3 = T^ 1 is the inverse temperature; A = e' 3 ^ denotes the fugacity and fj, the chemical potential; E is the particle 
energy equal to \J m 2 + p 2 with m being the particle mass and p its momentum; the upper sign is for fermions while 
the lower one for bosons. 

Since the inclusive azimuthal distribution is again given by (||), we get 




where 

^_ f d 3 p \~ 1 e /3E 

As seen, $ is an intensive thermodynamic quantity, i.e. it is independent of the system volume. It is also independent 
of the number of the particle internal degrees of freedom. One easily observes that $ < for fermions, $ > for 
bosons and $ = in the classical limit (A -1 S> 1). 
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When the particles are massless and their chemical potential vanishes (A = 1), $ can be calculated analytically and 
the result reads 



vnV 6 c(3)l 



/2/3 



1 



-0.082 
0.309 



where £(x) is the Riemann zeta function with £(3) = 1.202. 

In Fig. 1 we present with dashed lines the measure of 0— fluctuations in the ideal pion gas as a function of 
temperature. The pions are, of course, massive (m^ = 140 MeV). The calculations are performed for several values of 
the pion chemical potential. The chemical equilibrium corresponds to p — 0. As seen, $ grows with the temperature. 
In the classical limit, when p, — ► 0, the measure vanishes. 

It is a far going idealization to treat a fireball at freeze-out as an ideal gas of pions. A substantial fraction of the 
final state pions come from the hadron resonances. These pions do not 'feel' the Bosc-Einstein statistics at freeze-out 
and consequently the value of <& should be significantly reduced. We estimate the role of resonances as in our earlier 
papers |H^0). The spectrum of pions, which originate from the resonance decays, is not dramatically different than 
that given by the equilibrium distribution |2l}| . Therefore, we treat the fireball at freeze-out as a mixture of 'quantum' 
pions - those called 'direct' - and the 'classical' ones which come from the resonance decays. The $— measure is again 
given by Eq. (0) but the formulas flsjfiol) are modified as 



d 3 p 
(2^)3 



1 

X -l e 0E _ I 



+ \ r e 



-0E 



P 



d 3 p 
(2^)3 



— + \ r e-? E 



( A -l e /3£_ 1)2 



The parameter A r is chosen is such a way that the number of 'classical' pions equals the number of pions from the 
resonance decays. Thus, A,- is temperature dependent. In the actual calculations we have taken into account the 
lightest resonances: p(770) and w(782) which give the dominant contribution. The life time of p, which is 1.3 fm/c, 
is not much longer than the time of the fireball decoupling and some pions from the p decays can still 'feel' the 
effect of Bose statitistics. Therefore, the contribution of p to the 'classical' pions is presumably overestimated in our 
calculations. Since we neglect the heavier resonances and weakly decaying particles, which also contribute to the final 
state pions, the two effects partially compensate each other. In any case, our calculations show that the resonances 
do not change the value of $ dramatically in the domain of temperatures of interest. 

In Figs. 1 the solid lines represent $— measure which includes the resonances. The chemical potentials of p and uj 
are assumed to be equal to that of pions. As seen, the role of the resonances is negligible at the temperatures below 
100 MeV but above this temperature the resonances reduce the fluctuations noticeably. The freeze-out temperature 
in Pb-Pb collisions at 158 GeV per nucleon, which is obtained by means of the simultaneous analysis of the single 
particle spectra and the two-particle correlations, is about 120 MeV |22|. For T = 120 MeV and p = the <I>— measure 
equals 0.078, when the resonances are neglected, and is reduced to 0.066 when the resonances are taken into account. 
One observes that the effects of the quantum statistics and transverse flow are of comparable size. 

The $— measure given by Eq. (^) corresponds to the second moment of the fluctuating quantity. It has been 
suggested [^3| to use the higher moments in an analogous way. However, we have shown |2(| that only the third 
moment measure preserves the advantageous properties of <f> while the higher moment measures do not. We have also 
argued [^(J that the simultaneous usage of «3>2 and $3 may help in identifying the origin of correlations observed in the 
final state of heavy-ion collisions. Unfortunately, the third moment measure is useless in the studies of </>— fluctuations. 
One easily shows that due to the symmetry z 3 = and (Z 3 ) = when the variable x is identified with the azimuthal 
angle. 

We conclude our considerations as follows. The $— measure, which can be easily applied to the experimental 
data, seems to be a useful tool to analyze the azimuthal fluctuations in heavy-ion collisions. It is sensitive to the 
different nontrivial fluctuations, in particular those caused by the transverse flow and quantum statistics which have 
been quantitatively discussed here. The $— measure analysis combined with other techniques, such as the Fourier 
expansion method will help to study various sources of the fluctuations. 
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Figure Caption 



Fig. 1. measure of (^—fluctuations in the hadron gas as a function of temperature for four values of the chemical 
potential. The resonances are either neglected (dashed lines) or taken into account (solid lines). The most upper 
dashed and solid lines correspond to ^ = 70 MeV, the lower ones to fx = 0, etc. 
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